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ABSTRACT 

We present an analysis of the environmental dependence of bars and bulges in disc 
galaxies, using a volume-limited catalogue of 15810 galaxies at z < 0.06 from the 
Sloan Digital Sky Survey with visual morphologies from the Galaxy Zoo 2 project. 
We find that the likelihood of having a bar, or bulge, in disc galaxies increases when 
the galaxies have redder (optical) colours and larger stellar masses, and observe a 
transition in the bar and bulge likelihoods at = 2 x 10 10 Mq, such that massive 
disc galaxies are more likely to host bars and bulges. In addition, while some barred 
and most bulge-dominated galaxies are on the "red sequence" of the colour-magnitude 
diagram, we see a wider variety of colours for galaxies that host bars. We use galaxy 
clustering methods to demonstrate statistically significant environmental correlations 
of barred, and bulge- dominated, galaxies, from projected separations of 150kpc//i to 
3Mpc//i. These environmental correlations appear to be independent of each other: 
i.e., bulge-dominated disc galaxies exhibit a significant bar-environment correlation, 
and barred disc galaxies show a bulge-environment correlation. As a result of sparse 
sampling tests — our sample is nearly twenty times larger than those used previously — 
we argue that previous studies that did not detect a bar-environment correlation were 
likely inhibited by small number statistics. We demonstrate that approximately half 
of the bar-environment correlation can be explained by the fact that more massive 
dark matter haloes host redder disc galaxies, which are then more likely to have bars; 
this fraction is estimated to be 50 ± 10 % from a mock catalogue analysis and 60 ±5 % 
from the data. Likewise, we show that the environmental dependence of stellar mass 
can only explain a smaller fraction (25 ± 10 %) of the bar-environment correlation. 
Therefore, a significant fraction of our observed environmental dependence of barred 
galaxies is not due to colour or stellar mass dependences, and hence must be due to 
another galaxy property, such as gas content, or to environmental influences. Finally, 
by analyzing the projected clustering of barred and unbarred disc galaxies with halo 
occupation models, we argue that barred galaxies are in slightly higher-mass haloes 
than unbarred ones, and some of them (approximately 25%) are satellite galaxies in 
groups. We discuss the implications of our results on the effects of minor mergers and 
interactions on bar formation in disc galaxies. 

Key words: methods: statistical - galaxies: evolution - galaxies: structure - galaxies: 
spiral - galaxies: bulges - galaxies: haloes - galaxies: clustering - large scale structure 
of the universe 
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1 INTRODUCTION 

In recent years, there has been a resurgence in interest in 
the "secular" processes that could affect galaxy evolution 
(e.g., Weinzirl et al. 2009; Schawinski et al. 2010; Emsellem 
et al. 2011), driven by the growing understanding that ma- 
jor mergers are rare (e.g., Hopkins et al. 2010b; Darg et al. 
2010; Lotz et al. 2011), and may not play as important a role 
in galaxy evolution as had previously been thought (Parry 
et al. 2009; Dave et al. 2011). In particular, bars have been 
found to be common structures in disc galaxies, and are 
thought to affect the evolution of galaxies (e.g., Sellwood & 
Wilkinson 1993; Kormendy & Kennicutt 2004) and the dark 
matter haloes that host them (e.g., Debattista & Sellwood 
2000; Weinberg & Katz 2002). The abundance and proper- 
ties of barred galaxies have been analyzed in low and high- 
redshift surveys (e.g., Jogee et al. 2005; Sheth et al. 2005, 
2008; Barazza et al. 2008; Aguerri et al. 2009; Nair & Abra- 
ham 2010; Cameron et al. 2010; Masters et al. 2011; Hoyle 
et al. 2011; Ellison et al. 2011) and have been modeled with 
detailed numerical simulations, including their interactions 
with the host dark matter haloes (e.g., Valenzuela & Klypin 
2003; O'Neill & Dubinski 2003; Debattista et al. 2006; Heller 
et al. 2007; Weinberg & Katz 2007). 

Bars are extended linear structures in the central re- 
gions of galaxies, which form from disc instabilities and 
angular momentum redistribution within the disc (e.g., 
Athanassoula 2003; Berentzen et al. 2007; Foyle et al. 2008). 
Bars are efficient at driving gas inwards, perhaps sparking 
central star formation (e.g., Friedli et al. 1994; Ellison et al. 
2011), and thus may help to grow a central bulge component 
in galaxy discs (e.g., Dalcanton et al. 2004; Debattista et 
al. 2006; Gadotti 2011). Such bulges are sometimes referred 
to as "pseudo-bulges", to distinguish them from "classical" 
bulges, which are often thought to have formed from the 
hierarchical merging of smaller objects (e.g., Kormendy & 
Kennicutt 2004; Drory & Fisher 2007; De Lucia et al. 2011; 
Fontanot et al. 2011). 

Bars and (classical) bulges may also be related struc- 
tures and in some cases could form simultaneously. Galaxies 
with earlier-type morphologies, which have more prominent 
bulges, tend to have more, and longer, bars (Elmegreen & 
Elmegreen 1985; Weinzirl et al. 2009; Masters et al. 2011; 
Hoyle et al. 2011; Elmegreen et al. 2011; cf., Barazza et al. 
2008). In addition, at least in some galaxies, the bars and 
bulges have similar stellar populations (Sanchez-Blazquez 
et al. 2011). Nonetheless, there are some barred galaxies 
that lack bulges and many bulge-dominated galaxies that 
lack bars (e.g., Laurikainen et al. 2007; Perez & Sanchez- 
Blazquez 2011). 

Various classification methods have been developed to 
observationally identify bars, either visually or using auto- 
mated techniques, such as ellipse-fitting of isophotes and 
Fourier decomposition of surface brightness distributions 
(e.g., Erwin 2005; Aguerri et al. 2009; Gadotti 2009). These 
have yielded similar, but not always consistent, bar fractions 
(see discussions in Sheth et al. 2008; Nair & Abraham 2010; 
Masters et al. 2011). All bar identification methods are af- 
fected by issues such as inclination, spatial resolution, wave- 
length dependence, surface brightness limits, and selection 
biases (e.g., Menendez-Delmestre et al. 2007). 

In this paper, we use data from the Galaxy Zoo 2 project 



(see Masters et al. 2011), which provides detailed visual clas- 
sifications of ~ 250, 000 galaxies in the Sloan Digital Sky 
Survey (SDSS; York et al. 2000). Galaxy Zoo yields a rela- 
tively large catalogue of galaxies with reliable classifications 
in a variety of environments. It is particularly suited for 
analyses of the environmental dependence of the morpho- 
logical and structural properties of galaxies across a range 
of scales, as the large volume and sample size makes it less 
affected by cosmic variance than other catalogues. 

It has long been known that galaxy morphologies are 
correlated with the environment, such that spiral galaxies 
tend to be located in low-density regions and early-type 
galaxies in denser regions (e.g., Dressier 1980; Postman & 
Geller 1984; and confirmed by Galaxy Zoo: Bamford et al. 
2009; Skibba et al. 2009). There are a variety of ways to 
assess the correlation between galaxy properties and the en- 
vironment, such as fixed aperture counts and distances to 
nearest neighbors (see reviews by Haas et al. 2012; Muldrew 
et al. 2012). We follow Skibba & Sheth (2009) and Skibba 
et al. (2009) by using two-point galaxy clustering. 

There has been some recent work focused specifically 
on the environmental dependence of barred galaxies (van 
den Bergh 2002; Li et al. 2009; Aguerri et al. 2009; Mendez- 
Abreu et al. 2010; Giordano et al. 2011). All of these studies 
argue that there is little to no dependence of galaxy bars 
on the environment. Contrary to these results, Barazza et 
al. (2009) and Marinova et al. (2009, 2012) detect a slightly 
larger bar fraction in the cores of galaxy clusters, but of 
weak statistical significance, and Barway et al. (2011) find 
a higher bar fraction of faint SOs in group/cluster environ- 
ments. There is as yet no consensus on the environmental de- 
pendence of galaxy bars. These studies have been hampered 
by small number statistics, having typically a few hundred 
to a thousand galaxies at most. We improve upon this work 
by analyzing the environmental dependence galaxy bars and 
bulges in Galaxy Zoo 2, using a volume- limited catalogue 
consisting of 15810 disc galaxies in the SDSS. 

This paper is organized as follows. In the next section, 
we describe the Galaxy Zoo 2 data and our volume-limited 
catalogue. We introduce mark clustering statistics, and in 
particular, the marked correlation function, in Section [3] In 
Section [4] we show the distributions and correlations be- 
tween measures of bars and bulges. Then in Section [5] we 
present some of our main results, about the environmental 
dependence of barred and bulge-dominated galaxies, and we 
interpret the results with mock catalogues and halo occupa- 
tion models in Section [6] We end with a discussion of our 
results in Section [7] 



2 DATA 

2.1 Morphological Information from Galaxy Zoo 

To identify bars in local galaxies we use classifications pro- 
vided by members of the public through the Galaxy Zoo 
websitdM Specifically, we use classifications from the sec- 
ond phase of Galaxy Zoo (hereafter GZ2) which ran for 14 



http : //www . galaxyzoo . org 



© 0000 RAS, MNRAS 000, [Tlfil)] 



Environmental Dependence of Bars and Bulges 3 



months (between 9th Feb 2009 and 22nd April 20Wf\ In 
GZ2, volunteers were asked to provide detailed classifica- 
tions for the brightest (in terms of flux) 250,000 galaxies 
in the SDSS Data Release 7 (DR7; Abazajian et al. 2009) 
Main Galaxy Sample (MGS; Strauss et al. 2002). The se- 
lection criteria for GZ2 are m r < 17 and rgo > 3", where 
m r is the r-band Petrosian magnitude, and rgo is the radius 
containing 90% of this flux. Additionally, where the galaxy 
has a measured redshift, the selection 0.0005 < z < 0.25 is 
applied. 

Following the method of the original Galaxy Zoo (Lin- 
tott et al. 2008, 2011; hereafter GZ1), volunteers were asked 
to classify galaxies from the gri composite images. The com- 
plete GZ2 decision tree is shown in Fig. 1 of Masters et al. 
(2011). Users were presented each question in turn, with 
their progress down the tree depending on their previous 
answers. In the 14 months GZ2 ran, 60 million individual 
classifications were collected, with each galaxy in GZ2 hav- 
ing been classified by a median of 40 volunteers (i.e., 40 peo- 
ple answering the question at the top of the tree). As was 
discussed for GZ1 in Lintott et al. (2008) and Bamford et al. 
(2009), the conversion of these raw clicks into a unique classi- 
fication for each galaxy is a process similar to data reduction 
that must take into account possible spurious classifications 
and other problems. An iterative weighting scheme is used 
to remove the influence of unreliable users, and a cleaning 
procedure is applied to remove multiple classifications of the 
same galaxy by the same user. 

In what follows we will call the total number of cleaned 
and weighted classifications for a given question, Nx, and 
the fraction of positive answers to a given question px- For 
example, we will discuss N to tsd, which is the total (weighted) 
number of users classifying a given galaxy; features, the 
(weighted) fraction of such users identifying the galaxy as 
having features; A r b ar , the total (weighted) number of clas- 
sifications to the "bar question"; and pbar, the (weighted) 
fraction of users who indicated that they saw a bar. 

Examples of disc galaxies with a range of pbar are shown 
in Figure [l] The top row shows four randomly selected galax- 
ies with pbar = 0.0 (and a range of values of fracdeV, which 
is used to indicate the bulge size; see Sec. 2.2 for details). 
The lower rows show galaxies with larger values of pb ar ; the 
galaxies in the bottom row clearly have strong bars. 



2.2 Bulge Sizes 

While bulge size identification was present in the GZ2 clas- 
sification scheme (in the question of "How prominent is the 
central bulge?"), we choose instead to follow Masters et al. 
(2010a) and use the SDSS parameter "fracdeV" , which is 
a continuous indicator of bulge sizes in disc galaxies (see 
Kuehn & Ryden 2005; Bernardi et al. 2010) and is strongly 
correlated with the GZ2 bulge classification (Masters et al. 
2011). In the SDSS pipeline (Subbarao et al. 2002), galaxy 
light profiles are fit with both an exponential and de Vau- 
coulcurs profile (de Vaucouleurs 1948), and the model mag- 
nitude comes from the best fit linear combination of these 
two profiles. The parameter fracdeV indicates the fraction of 
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Figure 1. Example images of disc galaxies with a range of values 
of Pbar an <i fracdeV. The top row shows four randomly selected 
galaxies with Pbar = 0.0 (and a range of values of fracdeV, as 
indicated). The second row shows galaxies with Pb ar = 0.2, then 
Pbar = 0.5 in the row below, and finally Pbar = 1-0 in the bottom 
row. 



this model r-band magnitude that is contributed by the de 
Vaucouleur profile (Vincent & Ryden 2005). It is expected 
to have the value fracdeV= 1 in elliptical galaxies, and also 
bulge-dominated disc galaxies whose central light is domi- 
nated by a spheroidal bulge component. In pure disc galax- 
ies with no central light excess over an exponential disc, 
fracdeV= is expected. As we will see in Section [4] many 
galaxies have either fracdeV= or 1. 

The fracdeV parameter is likely to be most effective at 
identifying classical bulges, although any central excess of 
light over an exponential disc will drive fracdeV away from 
a zero value (Masters et al. 2010a). The Sersic index, which is 
closely related to fracdeV (Vincent & Ryden 2005), is corre- 
lated with the bulge-to-total luminosity ratio, but with some 
scatter (Gadotti 2009). Gadotti (2009) also shows that the 
Sersic index can be used to distinguish between classical and 
pseudo-bulges, although it cannot perfectly separate them 
(see also Graham 2011). The concentration rgo/rso, the ra- 
tio of the radii containing 90% and 50% of a galaxy's light in 
the r-band, is another morphological indicator (Strateva et 
al. 2001), and we have found that it exhibits a qualitatively 
similar clustering dependence as fracdeV. A comparison be- 
tween fracdeV, concentration, and GZ1 spiral and early- type 
classifications is shown in Masters et al. (2010a). 



2.3 Other Galaxy Properties 

In addition to morphological classifications and light profile 
shapes, we use various other parameters from the SDSS, in- 
cluding redshifts, (g — r) and (g — i) colours (from the model 
magnitudes), M r and Mi total magnitudes (for which we use 
the Petrosian magnitudes), and axial ratio, log(a/b) (from 
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the exponential model axial ratio fit in the r-band). All mag- 
nitudes and colours are corrected for Galactic extinction us- 
ing the maps of Schlegel, Finkbeiner & Davis (1998) and 
are A'-corrected to z = 0.0 using kcorrect v4_2 (Blanton 
& Roweis 2007). 

We compute stellar masses using the Zibetti et al. 
(2009) stellar mass calibration, which is based on the total 
magnitude °'°Mi and °'°(g — i) colour (the superscript "0.0" 
refers to the redshift of the ^-correction), with an absolute 
solar magnitude of M^q = 4.52 (Blanton et al. 2001), and 
assuming a Chabrier initial mass function (Chabrier et al. 
2003). We refer the reader to Zibetti et al. (2009) for details 
of the model. 

The stellar mass-to-light ratios typically have 0.2 — 
0.3 dex scatter. A more accurate method to estimate stellar 
masses would have been to apply stellar population mod- 
els (e.g., Maraston 2005) directly to the SDSS photometry 
in all five optical passbands. The Zibetti et al. (2009) cali- 
bration (which uses an updated version of Bruzual & Char- 
lot (2003) models) is nonetheless consistent with Maraston 
(2005), with no systematic offsets between their masses and 
with discrepancies only at young stellar ages, and is suffi- 
cient for the analysis in this paper. 



2.4 Volume-Limited Disc Galaxy Sample 

We perform our analysis on a volume-limited sample. Our 
catalogue is a subsample of the GZ2 catalogue, with limits 
-23.5 < °-°M r - 51og(/i) < -19.4 and 0.017 < z < 0.060. 
This catalogue is similar to the volume-limited catalogue 
used in Skibba et al. (2009; hereafter S09), but it has a 
slightly fainter absolute magnitude threshold because it is 
limited to slightly lower redshifts where we expect the bar 
identification in GZ2 to be most reliable (see Masters et 
al. 2011, hereafter Mil, and Hoyle et al. 2011 for further 
discussion of this choice); in addition, GZ2 has a slightly 
brighter flux limit than GZ1. 

The absolute magnitude threshold of our volume- 
limited catalogue approximately corresponds to M r < M* + 
1, where M* is the Schechter function break in the r- 
band luminosity function (Blanton et al. 2001). It also 
corresponds to an approximate stellar mass threshold of 
« 4 x 10 9 M Q (Zibetti et al. 2009) and a halo mass threshold 
of w 5 x 10 11 h' 1 M Q (Skibba & Sheth 2009), although there 
is substantial scatter between galaxy luminosity and stellar 
and halo masses. 

The absolute magnitude and redshift limits result in a 
catalogue of 45581 galaxies. We limit the sample further to 
log(a/6) < 0.3, which is approximately an inclination of 60°, 
in order to select face-on or nearly face-on galaxies. This is 
a comparable inclination cut to other recent studies of bars 
(e.g., Barazza et al. 2008; Sheth et al. 2008; Aguerri et al. 
2009) and identical to the cut used by Mil to study the bar 
fraction of GZ2. After this inclination cut, the catalogue is 
reduced to 32019 galaxies. 

We require a reasonable number of answers to the bar 
identification question in GZ2. As can been seen in the GZ2 
classification tree, in order to identify the presence of a bar 
in a galaxy, the volunteer must first identify the galaxy 
as "having features", and in addition answer "no" to the 
question "Could this be an edge-on disc?" We therefore 
limit the sample to TVbar > ^ tota i/4 (which is equivalent 
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Figure 2. Distribution in redshift and Petrosian r-band absolute 
magnitude, with the selection criteria of the volume-limited cat- 
alogue: -23.5 < °-°M r - 51og(/i) < -19.4 and 0.017 < z < 0.060. 
Without the N^,^ classification cut (selecting galaxies with suffi- 
cient classifications) and a/b axis ratio cut (selecting nearly face- 
on disc galaxies), the magnitude-redshift distribution is virtually 
identical; these cuts do not bias the catalogue. 



to PfoaturcsPnotcdgc-on > 0.25), resulting in a catalogue of 
15989 galaxies. We also remove a small number (179) of ob- 
jects with iVbar < 10 which may have bar identifications 
dominated by a small number of classifiers. 

Our resulting volume-limited catalogue comprises 15810 
nearly face-on disc galaxies with reliable bar classifications. 
The galaxy distribution in redshift and magnitude, and the 
cuts used to define the catalogue, are shown in Figure [2] 
Note that, because of the A r b ar cuts, elliptical galaxies are 
excluded from the sample. We emphasize that only disc 
galaxies (i.e., spiral galaxies and SOs) constitute the sample, 
including "bulge-dominated" disc galaxies and "pure disc" 
galaxies, to which we often refer as "disc-dominated" galax- 
ies. We will distinguish between bulge-dominated and disc- 
dominated galaxies with the fracdeV parameter (described 
in Sec.[2~2|. 



In principle, some of the selection criteria, such as the 
NhnT requirements, could bias our results by excluding cer- 
tain galaxies in an environmentally dependent way. Never- 
theless, Mil have tested this by comparing the luminos- 
ity, colour, axial ratio, and redshift distributions with and 
without this requirement, and have found no significant dif- 
ferences. Therefore, it is unlikely that there are any signifi- 
cant biases introduced; on the contrary, these criteria should 
eliminate biases by removing contaminating objects, such as 
elliptical galaxies or mergers with unreliable bar classifica- 
tions. 

We have also tested our clustering measurements as a 
function of the axial ratio a/b, and confirmed that they are 
not affected by the inclination cut. In particular, the corre- 
lation functions for different inclinations are within 0.03 dex 
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(7%), well within the error bars, and the mark correlation 
functions, described in the next section, are within 2%, ex- 
cept at separations of r p < 500 kpc//i, where they still agree 
within 10%. 

In addition, "fiber-collided" galaxies are not included 
in the catalogue. The thickness of the spectroscopic fibers 
means that some galaxies closer than 55" on the sky will 
be missing spectra and redshifts. This fiber-collision con- 
straint is partly alleviated by the fact that neighbouring 
plates have overlap regions, but it still results in 7% of tar- 
geted galaxies not having a measured redshift (Zehavi et al. 
2005) and could significantly affect clustering measurements, 
especially at separations smaller than 100 kpc/h (Guo, Ze- 
havi & Zheng, 2011). Nevertheless, we focus our analysis 
on marked correlation functions, in which the effects of fiber 
collisions are expected to cancel out (see Eqn.[5]in Section|3j 
where we describe the marked correlation functions). More- 
over, the fiber assignments were based solely on target po- 
sitions, and in cases where multiple targets could only have 
a single fiber assigned, the target selected to be observed 
was chosen randomly — hence independently of galaxy prop- 
erties. Therefore, we argue that the effects of fiber collisions 
are likely to be negligible for the marked correlation func- 
tions. 

Throughout this paper we assume a spatially flat cos- 
mology with Q m — 0.3 and J1a = 1 — fim- We write the 
Hubble constant as Ho = WOh km s _1 Mpc~ . 



3 MARK CLUSTERING STATISTICS AND 
ENVIRONMENTAL CORRELATIONS 

We characterize galaxies by their properties, or "marks", 
such as their luminosity, colour, morphological type, stel- 
lar mass, star formation rate, etc. In most galaxy clustering 
analyses, a galaxy catalogue is cut into subsamples based 
on the mark, and the two-point clustering in the subsample 
is studied by treating each galaxy in it equally (e.g., Madg- 
wick et al. 2003, Zehavi et al. 2005, Tinker et al. 2008). 
These studies have shown that galaxy properties are corre- 
lated with the environment, such that elliptical, luminous, 
and redder galaxies tend to be more strongly clustered than 
spiral, fainter, and bluer galaxies. 

Nonetheless, the galaxy marks in these studies are used 
to define the subsamples for the analyses, but are not consid- 
ered further. This procedure is not ideal because the choice 
of critical threshold for dividing galaxy catalogues is some- 
what arbitrary, and because throwing away the actual value 
of the mark represents a loss of information. In the current 
era of large galaxy surveys, one can now measure not only 
galaxy clustering as a function of their properties, but the 
spatial correlations of the galaxy properties themselves. We 
do this with "marked statistics" , in which we weight each 
galaxy by a particular mark, rather than simply count galax- 
ies as "one" or "zero". 

Marked clustering statistics have been applied to a vari- 
ety of astrophysical datasets by Beisbart & Kerscher (2000), 
Gottlober et al. (2002), and Martinez et al. (2010). Marked 
statistics are well-suited for identifying and quantifying cor- 
relations between galaxy properties and their environments 
(Sheth, Connolly & Skibba 2005). They relate traditional 
unmarked galaxy clustering to the clustering in which each 



galaxy is weighted by a particular property. Marked statis- 
tics are straightforward to measure and interpret: if the 
weighted and unweighted clustering are significantly differ- 
ent at a particular scale, then the galaxy mark is correlated 
(or anti-correlated) with the environment at that scale, and 
the degree to which they are different quantifies the strength 
of the correlation. In addition, issues that plague traditional 
clustering measurements, such as incompleteness and com- 
plicated survey geometry, do not significantly affect mea- 
surements of marked statistics, as these effects cancel out to 
some extent, since the weighted and unweighted measure- 
ments are usually similarly affected. Mark correlations have 
recently been measured and analyzed in galaxy and dark 
matter halo catalogues (e.g., Sheth & Tormen 2004; Sheth 
et al. 2006; Harker et al. 2006; Wetzel et al. 2007; Mateus 
et al. 2008; White & Padmanabhan 2009; S09). Finally, the 
halo model framework has been used to interpret the cor- 
relations of luminosity and colour marks in terms of the 
correlation between halo mass and environment (Skibba et 
al. 2006, Skibba & Sheth 2009). We focus on morphological 
marks here, in particular, the likelihood of galaxies having 
a bar or bulge component. 

There are a variety of marked statistics, but the easiest 
to measure and interpret is the marked two-point correlation 
function. The marked correlation function is defined as the 
following: 



M(r) 



W(r) 



(1) 



where £(r) is the two-point correlation function, the sum 
over galaxy pairs separated by r, in which all galaxies are 
"weighted" by unity. W(r) is the same sum over galaxy pairs 
separated by r, but now each member of the pair is weighted 
by the ratio of its mark to the mean mark of all the galax- 
ies in the catalogue (e.g., Stoyan & Stoyan 1994). That is, 
for a given separation r, £(r) receives a count of 1 for each 
galaxy pair, and W(r) receives a count of WiWj for W(r). 
The fact that the real-space (not redshift-distorted) marked 
statistic M(r) can be approximately estimated by the sim- 
ple pair count ratio WW/ DD (where DD are the counts of 
data-data pairs and WW are the weighted counts), with- 
out requiring a random galaxy catalogue, implies that the 
marked correlation function is less sensitive than the un- 
marked correlation function to the effects of the survey edges 
(Sheth, Connolly & Skibba 2005). In effect, the denominator 
in Eqn. [T] divides out the contribution to the weighted cor- 
relation function which comes from the spatial contribution 
of the points, leaving only the contribution from the fluctua- 
tions of the marks. The mark correlation function measures 
the clustering of the marks themselves, in environments of 
a given scale. 

In practice, in order to obviate issues involving redshift 
distortions, we use the projected two-point correlation func- 
tion 



w p (r p ) = 



[ dr£(r p ,7r) = 2 f 

J J r„ 



dr ■ 



r£(r) 



(2) 



where r = -y/r p 2 + n 2 , r p and -k are the galaxy separations 
perpendicular and parallel to the line of sight, and we inte- 
grate up to line-of-sight separations of ir = 40Mpc//i. We 
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estimate £(r p , 71") using the Landy & Szalay (1993) estimator 
DD - 2DR + RR 



RR 



(3) 



where DD, DR, and RR are the normalized counts of data- 
data, data-random, and random-random pairs at each sep- 
aration bin. Similarly, the weighted projected correlation 
function is measured by integrating along the line-of-sight 
the analogous weighted statistic 

WW - 2WR + RR 



W(r p ,ir) = 



RR 



(4) 



where W refers to a galaxy weighted by some property (for 
example, pb ar or fracdeV; see Section [4b, and R now refers 
to an object in the catalogue of random points, weighted by 
a mark chosen randomly from its distribution. 

We then define the marked projected correlation func- 
tion: 



M p (r p ) 



1 + W p (r p )/r p 
1 + w p (r p )/r p 



(5) 



which makes M p (r p ) w M(r) on scales larger than a few 
Mpc, in the linear regime. The projected correlation func- 
tions w p and W p (r p ) are normalized by r p , so as to be made 
unitless. On large scales both the real-space and projected 
marked correlation functions (Eqns [I] and [5| will approach 
unity, because at increasing scale the correlation functions 
£(r) and W(r) (or w p (r p ) and W p (r p )) become small as 
the universe appears nearly homogeneous. The simple ra- 
tio of the weighted to the unweighted correlation function 
W(r)/£(r) (or W p (r p ) /w p (r p )) approaches unity similarly, 
provided that there are sufficient number statistics and the 
catalogue's volume is sufficiently large. 

For the correlation functions and error measurements, 
which require random catalogues, we use the hierarchi- 
cal pixelization scheme SDSSPijj^] which characterizes the 
survey geometry, including edges and holes from miss- 
ing fields and areas near bright stars. This pixelization 
scheme has been used for clustering analyses (Scranton 
et al. 2005, Hansen et al. 2009) and lensing analyses 
(Sheldon et al. 2009). We use the Scranton et al. code, 
j ack_random_polygon, to construct the catalogues, and we 
use at least twenty times as many random points as in the 
data for all of the clustering measurements. 

We estimate statistical errors on our measurements us- 
ing "jack-knife" resampling. We define iV S ub = 30 spatially 
contiguous subsamples of the full dataset, and the jack-knife 
subsamples are then created by omitting each of these sub- 
samples in turn. The scatter between the clustering mea- 
surements from the jack-knife samples is used to estimate 
the error on the clustering statistics, w p , W p , and M p . The 
jack-knife covariance matrix is then, 



Covai(xi, Xj 



N sub 



N sul 



« 3Ub 



%i ) yEj 



(6) 



k=l 



where Xi is the mean value of the statistic x measured in the 
i radial bin in all of the samples (see Zehavi et al. 2005; 
Norberg et al. 2009). As shown by these authors (see also 
McBride et al. 2011), however, the jack-knife technique only 
recovers a noisy realization of the error covariance matrix, 



http: //dls .physics .ucdavis . edu/'scranton/SDSSPix 



as measured from mock catalogues, but in any case, our 
results are not sensitive to correlated errors in the clustering 
measurements. 



4 RESULTS: DISTRIBUTIONS AND 

CORRELATIONS OF BAR AND BULGE 
PROPERTIES 

The structural galaxy properties that we examine in this 
paper are the bar fraction or probability, pbar, and fracdeV, 
which quantifies bulge strength (see Sections 2.1 and 2.2|. 
Using Pbar, we can compute the bar fraction of galaxies in 



the volume-limited catalogue, such as with a pbar threshold, 
though as noted by Mil, Galaxy Zoo tends to identify bars 
that are consistent with optically-identified strong bars (i.e., 
SB types). If one generously counts barred galaxies as those 
with pbar > 0.2, the bar fraction is /bar = 48.8±0.5% (where 
the error is estimated with bootstrap resampling), while if 
one counts those with pbar > 0.5, the fraction is 25.3 ± 
0.4%. The latter can be compared to Mil, who obtain a 
fraction of 29.4 ± 0.5% for their sample. Our slightly lower 
fraction may be due to the different selection criteria, such 
as the fact that the Mil sample has slightly more lower 
luminosity galaxies. For a comparison of bar classifications 
in GZ2, RC3, Barazza et al. (2008), and Nair & Abraham 
(2010), we refer the reader to Masters et al. (2012). 

In this section and Section [5] we will analyze pb ar and 
fracdeV, their environmental dependence, and their relation 
to each other and to galaxy colour and stellar mass. We will 



later (Sec. 5.2 and 6.3 1 compare the clustering of barred and 



unbarred disc galaxies. At that point, we separate the barred 
and unbarred galaxies as those with pbar > 0.2 and < 0.2, 
respectively, because it approximately splits the sample in 
half. The pbar > 0.5 threshold appears to identify strong 
bars, while we associate the range 0.2 < pbar < 0.5 with 
weak bars (Fig. [TJ . 

We first show the pbar and fracdeV distributions Fig- 
ure [3] The pbar distribution is smooth, with most galax- 
ies in the catalogue (62%) having pbar < 0.3. In contrast, 
the fracdeV distribution is peaked near and 1; namely, 
most galaxies are either distinctly disc-dominated or bulge- 
dominated. Recall though that, because of our selection cri- 
teria, all of the bulge-dominated galaxies in the catalogue 
have spiral arms or discs (i.e., elliptical galaxies are ex- 
cluded). We have tested that the clustering dependence of 
fracdeV is not very sensitive to its distribution: rescaling it 
to have a smoother distribution (such as that of pbar) yields 
clustering measurements within 5% of the results shown in 
the next section. 

As discussed in Section [l] some authors have argued 
that the formation and evolution of bars and bulges could 
be related, depending on the type of bulge, gas content, 
and angular momentum distribution (e.g., Debattista et al. 
2006; Laurikainen et al. 2007). Nevertheless, we find that 
Pbar and fracdeV are not simply, or monotonically, corre- 
lated, as can be seen from the distribution of pbar versus 
fracdeV in Figure [4] We find that a large fraction of bulge- 
dominated galaxies are barred (in the upper right corner 
of the figure) and a large fraction are unbarred (upper left 
corner): for example, of those with fracdeV> 0.7, 27% have 
Pbar > 0.7 and 49% have Pbar < 0.3. On the other hand, 
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Figure 3. Distribution of bar likelihood (pban blue dashed his- 
togram) and bulge strength (fracdeV, red solid histogram) of 
galaxies in the volume-limited catalogue. The histograms have 
been slightly offset from each other, for clarity. 




0.2 0.4 0.6 0.8 

Pbar 



Figure 4. Distribution of Pbar versus fracdeV, with contours in- 
dicating 10, 20, 33, 50, and 75% of the maximum counts. 

disc-dominated galaxies are mostly unbarred (76% of those 
with fracdeV< 0.3 have pbar < 0.3), and only a few per cent 
have bars — the lower right quadrant of the figure is empty. 
These results are consistent with Mil, who showed that the 
bar fraction of disc galaxies increases with fracdeV, which 
is clearly the case for galaxies with pbar > 0.5 on the right 
half of the figure. 

We show the colour-magnitude distribution of the cata- 



Figure 5. The optical colour-magnitude diagram (CMD), using 
extinction-corrected g—r colour and r-band magnitude, with filled 
gray contours indicating 5, 10, 20, 50, and 75% of the maximum 
counts. The majority of the barred disc galaxies (withpbar > 0-5) 
and bulge-dominated disc galaxies (with fracdeV> 0.9) are in- 
dicated by the blue and red contours, respectively. Pbar > 0.5 
and fracdeV> 0.9 select a similar number of galaxies in the 
catalogue. The black line indicates the red sequence separator 
(g - r = 0.63 - 0.02 (M r - 51og(/i) + 20)), used to identify red 
spirals in Masters et al. (2010b). 



logue in Figure [5] Many of the galaxies in the catalogue are 
disc-dominated, and the majority of them are located in the 
"blue cloud" , the bluer mode of the bimodal colour distri- 
bution (e.g., Skibba 2009). Applying the colour-magnitude 
separator used for red spiral galaxies (Masters et al. 2010b; 
similar to that of S09), we find that only 31% of the galax- 
ies in the whole sample are on the red sequence, while 74% 
of the bulge-dominated galaxies (fracdeV> 0.9, red contour) 
are on the red sequence (and are not highly inclined, so they 
are likely to have older stellar populations, rather than being 
dust reddened). Barred galaxies (pbar > 0.5, blue contour) 
have a much wider range of colours, with 48% of them on 
the red sequence. In addition, barred galaxies have a bimodal 
colour distribution based on bar length, such that galaxies 
with longer bars are on the red sequence (Hoyle et al. 2011). 

Next, in Figure[6] we show the median (and 1 1% and 1-er 
ranges) of pbar and fracdeV as a function of colour and stellar 
mass. The majority of blue galaxies and low-mass galaxies 
are disc-dominated, while most red and massive galaxies are 
bulge-dominated. The bar probability is also positively cor- 
related with colour and stellar mass, consistent with other 
studies (e.g., Sheth et al. 2008; Nair & Abraham 2010; Mil), 
such that redder and more massive galaxies are more likely 
to have bars. Some studies have found a bimodal distribu- 
tion of bars, such that blue, low-mass, or Sc/Sd-type galaxies 
also often have bars (Barazza et al. 2008; Nair & Abraham 
2010), in contrast to our finding that the median pb ar < 0.2 
at bluer colours and lower stellar masses. It is possible that 
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Figure 6. Left: Pbar (blue dashed line) and fracdeV (red solid line) as a function of g — r extinction-corrected colour. Right: Pbar (blue 
dashed line) and fracdeV (red solid line) as a function of stellar mass. Running medians are shown as the thicker lines, the 1 /2-a range is 
shown by the thinner lines, and the hatched regions indicate the 1-cr range between the 16 and 84 percentiles. As discussed in the text, 
Pbar an d fracdeV appear to transition at a similar colour and mass scale: g — r 0.6 and M* S3 2 X 10 10 Mq. 



in Galaxy Zoo, a large fraction of weak or short bars are 
missed in these galaxies. 

Compared to the correlation with fracdeV, the corre- 
lations with pbar in Figure [6] are not as strong and have 
more scatter, especially at the red and massive end. In other 
words, red and massive galaxies are more likely to have bars 
than blue and less-massive galaxies, but nonetheless there 
are many red and massive galaxies that lack bars. Either 
these galaxies never formed bars, or perhaps more likely, it 
is possible that they had bars in the past that were weakened 
(so that they are no longer detectable by GZ2) or destroyed; 
some galaxies may even have multiple episodes of bar for- 
mation in their lifetime (Bournaud & Combes 2002). 

It is interesting that the transition from mostly un- 
barred to mostly barred galaxies and from disc-dominated 
to bulge-dominated galaxies occurs at similar colours and 
stellar masses. The colour transition occurs at extinction- 
corrected g — r ~ 0.6, in the "green valley" of the colour- 
magnitude distribution (e.g., Wyder et al. 2007), between 
the blue and red peaks of the distribution (see Figure pi). 
The stellar mass transition occurs at M» ~ 2 x 1O 1O M0, 
and is similar to the mass scale identified by Kauffmann et 
al. (2003; see also Schiminovich et al. 2007), above which 
galaxies have high stellar mass surface densities, high con- 
centration indices typical of bulges, old stellar populations, 
and low star formation rates and gas masses. 



5 RESULTS: CLUSTERING OF GALAXIES 
WITH BARS AND BULGES 

We now explore the environmental dependence of disc galax- 
ies with bars and bulges by measuring marked projected cor- 



relation functions (described in Section j3j). For the marks, 
we use pbar and fracdeV, which indicate the presence or lack 
of a bar and bulge, respectively. We first present the total 
environmental dependence of these two galaxy properties in 
Section [5. 1| and then we attempt to separate their environ- 
mental dependences in Section [5^2] 

5.1 Environmental Correlations of Bar and Bulge 
Probability 

In the upper panels of Figure [7] we show the projected clus- 
tering of the galaxies in our (disc-dominated) catalogue. In 
the lower panels, we show the the marked correlation func- 
tions, quantifying the environmental correlations of pbar and 
fracdeV across a wide range of scales. The errors of the mea- 
surements are estimated using jack-knife resampling, and are 
analyzed in more detail in Appendix [A] As noted in the ap- 
pendix, a significant fraction of the error estimates is due to 
a single outlying jack-knife subsample. 

We see a statistically significant environmental correla- 
tion in Figure [7] for both pbar and fracdeV, which means 
that disc galaxies with bars and those with bulges tend 
to reside in denser environments on average. The environ- 
mental correlation is especially strong at small scales, at 
r p < 2Mpc//i. At these scales, the clustering signal is domi- 
nated by the "one-halo term" (pairs of galaxies within dark 
matter haloes), while at larger spatial scales the "two-halo 
term" (pairs of galaxies in separate haloes) dominates the 
clustering (e.g., Zehavi et al. 2004). In Section [5] we will 
interpret these correlation functions with the halo model of 
galaxy clustering (see Cooray & Sheth 2002 for a review). 

It is important to quantify the statistical significance 
of the M(r p ) measurements, the degree to which they are 
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Figure 7. Upper panel: projected correlation function w p (r p ) (circle points) and weighted projected correlation functions W p (r p ). Lower 
panel: marked projected correlation functions, using Pbar niark (left) and fracdeV mark (right). Recall that the marked correlation 
function is denned as M(r p ) = [1 + W p (r p ) /r p ]/[l + w p (r p )/r p ] (Eqn.[5|. 



inconsistent with unity. (A result of unity occurs when the 
weighted and unweighted correlation functions are the same, 
i.e., when the weight is not correlated with the environ- 
ment.) Since the errors are correlated, the statistical sig- 
nificance should be quantified using the covariance matrices 
(see Eqn. [6] and Appendix [A]): 

alt - (M - 1) T Covar 1 (M - 1), (7) 

where M is the pbar or fracdeV mark correlation function, 
and M — 1 is its deviation from unity. (This is similar to 
the way one would compute the x 2 of a fit, where in this 
case a good "fit" with low af.j would be a measurement 
consistent with unity, or no environmental correlation). The 
Pbar and fracdeV mark correlation functions yield a\ t = 39.9 
and 43.4, respectively, which could be interpreted as 6.3 and 
6.6a significance for the marks. Note that this significance 
estimate is not highly dependent on the binning: narrower 
r p bins, for example, would yield more measurements with 
M(r p ) > 1, but they would have larger and more correlated 
errors. 

The positive correlation at particular spatial separa- 
tions r p implies that galaxies with larger values o/pbar and 
fracdeV tend to be located in denser environments at these 
scales. This is one of the main results of the paper. This 
can also be seen in the upper panels of the figures, in which 
the weighted correlation functions are larger than the un- 
weighted ones. 

The environmental dependence of bulges is not surpris- 
ing; because of the "morphology-density relation" , in which 
dominant bulge components are associated with earlier-type 
morphologies, it is expected that bulge-dominated galaxies 
tend to be located in denser environments, in groups and 
clusters (e.g., Postman & Geller 1984; Bamford et al. 2009). 
Similar two-point clustering analysis have also clearly shown 



that bulge-dominated disc galaxies tend to be more strongly 
clustered than disc-dominated ones, on scales of up to a few 
Mpc (Croft et al. 2009; S09). Semi-analytic models, using 
stellar bulge-to-total ratios, have similarly predicted that 
bulge-dominated disc galaxies tend to form in more mas- 
sive dark matter haloes (e.g., Baugh et al. 1996; Benson & 
Devereux 2010; De Lucia et al. 2011). Many theorists have 
argued that bulge formation is linked to minor and major 
galaxy mergers, and mergers and interactions tend to be 
more common in denser environments, especially in galaxy 
groups (e.g., Hopkins et al. 2009; Hopkins et al. 2010a; Mar- 
tig et al. 2012); there is some observational evidence in favor 
of the link between bulges and mergers as well (e.g., Ellison 
et al. 2010). 

On the other hand, one might not expect barred galax- 
ies to be correlated with the environment, if bars form en- 
tirely by internal secular processes. Some recent studies have 
argued that barred and unbarred galaxies are located in sim- 
ilar environments, or have only weak evidence that barred 
galaxies are more strongly clustered at small scales (Mari- 
nova et al. 2009; Li et al. 2009; Barazza et al. 2009; Martinez 
& Muriel 2011; Wilman & Erwin 2012). Nevertheless, our 
volume-limited GZ2 catalogue is much larger than those of 
these studies (which except for Martinez & Muriel consist 
of ~ 1000 galaxies, or ~ l/17 th as many as ours), and as 
stated in Section [3j an advantage of mark clustering statis- 
tics is that one can analyze the entire catalogue, without 
splitting it and without requiring a classification of "clus- 
ter", "group", and "field" environments. We are thus able 
to quantify the correlation between bars and the environ- 
ment with greater statistical significance. 

To test the effect of small number statistics, we per- 
formed a number of sparse sampling measurements (ran- 
domly selecting galaxies, or only selecting galaxies in small 
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subregions or redshift slices of the sample). In general, we 
find that if fewer than ~ 2000 galaxies were used for our clus- 
tering measurements, we would not have a statistically sig- 
nificant detection of the environmental dependence of pbar, 
and the unmarked correlation functions would be too noisy. 
This may explain why previous studies of smaller galaxy 
catalogs did not detect a significant bar-environment corre- 
lation. 

The pbar marked correlation function in Figure [7J 
Mpbar(T" p ), increases in strength with decreasing spatial sep- 
aration. Such a trend is expected when the mark is positively 
correlated with the environment. An exception to this corre- 
lation occurs at small separations (r p ~ 100 kpc/h), where 
the result is consistent with no correlation at all. Fiber colli- 
sions could affect the lack of correlation at these scales, but 
at most 7% of the targeted galaxies lack measured redshifts, 
and we obtain M Pbar (100 kpc/h) ~ 1 for all jack-knife sub- 
samples, so this low mark correlation is likely a real effect. 
The weakening environmental correlation at small separa- 
tions suggests that whatever conditions that make bars more 
likely at larger scales are removed in close pairs of galax- 
ies. For example, close pairs are more likely to experience 
mergers, and bars may be weakened or destroyed immedi- 
ately following merger activity, although new bars may form 
later (Romano-Diaz et al. 2008). Nonetheless, Marinova et 
al. (2009) and Barazza et al. (2009) find that bar fractions 
are slightly larger in cluster cores, although this is of weak 
statistical significance according to the authors. More re- 
cently, Nair & Ellison (in prep.) find that the bar fraction 
of disc galaxies decreases as pair separation decreases, con- 
sistent with our results. 

It is interesting that Mpbar()>) peaks at approximately 
400 kpc/h (more precisely, the bin's range is 316 < r p < 
525 kpc/h). Many of the galaxies contributing to the signal 
at these scales are likely "satellite" galaxies in groups, rather 
than the central galaxies. In fact, considering that these are 
disc galaxies and that pbar is correlated with colour, it is 
likely that many of these are the same objects as the "red 
spirals" discussed in S09 (most of which have bars, according 
to Masters et al. 2010b; Mil), a relatively large fraction of 
which are satellites (/ sa t ~ 1/3). 

We also show (unmarked) clustering of barred versus 
non-barred galaxies (pbar > 0.2 and < 0.2), and bulge- 
dominated versus disc-dominated ones (fracdeV> 0.5 and 
< 0.5), in the upper panels of Figure [8] At large scales 
(r v > 2Mpc//i), their clustering strength is the same. 
At smaller separations, however, barred galaxies tend to 
be more strongly clustered than unbarred ones and bulge- 
dominated galaxies tend to be more strongly clustered than 
disc-dominated ones. The scale at which the correlation 
functions diverge corresponds to the scale of the transition 
from the "one-halo term" (pairs of galaxies within haloes) 
to the "two- halo term" (galaxies in separate haloes). These 
clustering measurements then suggest that barred and un- 
barred galaxies may reside in the same dark matter haloes, 
but the former are more likely to be central galaxies than the 
latter. The same applies for the presence/absence of bulges 
in central/satellite galaxies. We will return to this issue when 
we apply halo occupation modeling to the measured corre- 
lation functions, in Section [6. 3| 



5.2 Disentangling the Environmental Correlations 

As we have shown in previous sections, disc galaxies with 
large bulges are more likely to have a bar (see Fig. Hk. We 
have also shown that both bulge-dominated discs and discs 
with bars are more strongly clustered than average (Fig. [7b. 
We address in this Section the question of whether one of 
these two galaxy properties is more dependent on the envi- 
ronment, or whether their environmental correlations are in- 
dependent. That is to say, we will determine whether bulge- 
dominated galaxies with bars are more strongly clustered 
than bulge-dominated galaxies without bars, and whether 
barred galaxies with bulges are more strongly clustered than 
barred galaxies with no or small bulges. 

In addition, we know that disc galaxies hosting bars 
tend to be redder and have higher stellar masses than those 
with weak or no bars (Fig. |6| . We will later address in Sec- 
tion |6.1| whether the environmental correlations of galaxy 
colour or stellar mass (e.g., Skibba & Sheth 2009; Li & White 
2009) can account for the environmental correlation we have 
observed of bars. 

In the lower panels of Figure [8] we show the pb ar 
mark correlation functions for bulge-dominated and disc- 
dominated galaxies (fracdeV> 0.5 and < 0.5), as well as the 
fracdeV mark correlation functions for barred and unbarred 
galaxies (pbar > 0.2 and < 0.2). Using the fracdeV> 0.5 
threshold, 44% of our (disc) galaxy catalogue is bulge- 
dominated, and using pb ar > 0.2, 49% of it is barred. Fol- 
lowing the procedure described in the appendix of S09, the 
mark correlations are shown when the marks are rescaled so 
that they have the same distribution as that of the whole 
sample (see Fig.|3|. Such a rescaling is necessary in order to 
compare the mark correlation functions. (In this case, the 
mark correlation measurements are similar, within ~10%, 
when the mark distributions are not rescaled.) 

The pbar and fracdeV mark correlation functions are all 
still above unity, but they are statistically significant only 
for bulge-dominated (fracdeV> 0.5) and barred (pb ar > 0.2) 
galaxies, respectively. Using Eqn. [7J these pb ar and fracdeV 
mark correlations both have a statistical significance of 
6<t. In other words, bulge-dominated galaxies exhibit a sig- 
nificant bar-environment correlation, and barred galaxies 
exhibit a bulge-environment correlation. Considering that 
these residual environmental correlations are so significant, 
it appears that the environmental dependencies of barred 
and bulge-dominated galaxies are somewhat independent of 
each other: the bar-environment correlation is not due to 
the bulge-environment correlation, and vice versa. (The en- 
vironmental dependencies of bars and pseudo-bulges, how- 
ever, may be more closely related, as discussed in the in- 
troduction.) Lastly, we point out that though SOs are to 
some extent environmentally dependent (Hoyle et al. 2011; 
Wilman & Erwin 2012), they are not likely to be driving the 
bar-environment correlation in the left panels of Fig. [8] SOs 
do not have a particularly large bar fraction compared to 
their spiral counterparts (Laurikainen et al. 2009; Buta et 
al. 2010), and their bar fraction does not exhibit a significant 
environmental dependence (Barway et al. 2011; Marinova et 
al. 2012). 
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Figure 8. Left plot, lower panel: Pbar mark correlation functions with fracdeV< 0.5 (disc-dominated; open triangle points) and > 0.5 
(bulge-dominated; solid triangles). Right plot, lower panel: fracdeV mark correlation functions with Pbar < 0.2 (unbarred; open triangles) 
and > 0.2 (barred; solid triangles); the split is done at 0.2 in order to have a similar number of galaxies in each subsample. Mark 
distributions have been rescaled to match distribution of whole sample (Fig. [3j . For clarity, the points are slightly offset from each other, 
and the points are connected for bulge-dominated and barred galaxies, to guide the eye. Upper panels show the (unweighted) correlation 
functions, for the full sample (same as in Fig.[7]l and the subsamples, indicated by the legends in the figures. 



6 INTERPRETATION OF THE 

Pbar-ENVIRONMENT CORRELATION 



In the previous section, we quantified the environmental de- 
pendence of barred galaxies, using projected clustering mea- 
surements with the largest catalogue of galaxies with bar 
classifications to date. Here we perform tests and analyses 
of these results, in order to better understand the origin of 
these environmental correlations. 

We also quantified the environmental dependence of 
galaxy bulges, but as stated in Section |5.1| this has been 
thoroughly studied already and is closely related to the 
morphology-density relation. Furthermore, the colour and 
stellar mass dependence of the morphology-density relation 
has been studied elsewhere (e.g., Kauffmann et al. 2004; 
Blanton et al. 2005; Park et al. 2007; van der Wei et al. 
2010), including with Galaxy Zoo data (Bamford et al. 2009; 
Skibba et al. 2009), so we will not study it further here. 

In Section [6. 1[ we examine the stellar mass and colour 
dependence of the measured pbar-environment correlation. 
Then in Section [6. 2[ we use mock galaxy catalogues to pre- 
dict the strength of the pbar-environment correlation if it 
were entirely due to redder galaxies occupying more mas- 
sive dark matter haloes. Lastly, we apply halo occupation 
models to the clustering of barred and unbarred galaxies in 
Section 16.31 



6.1 Dependence of the environmental correlation 
on stellar mass and colour 

The probability of a galaxy being barred is correlated with 
its stellar mass (see Fig. [6)3; Nair & Abraham 2010), so it is 
important to ask whether the environmental dependence of 
barred galaxies measured in Section [5.1| is due to the envi- 
ronmental dependence of stellar mass. Li et al. (2009) argue 
that in their catalogue, at fixed stellar mass, the projected 
clustering of barred and unbarred galaxies is similar. With 
a much larger volume-limited catalogue, we can now ana- 
lyze the stellar mass and colour dependence with greater 
accuracy. 



6.1.1 Mark shuffling test 

Rather than splitting our catalogue by stellar mass and then 
measuring the pbar mark correlation function of the subcata- 
logues, we can take better advantage of the number statistics 
with a different test. Our procedure is as follows. We ran- 
domly shuffle the pbar marks at a given stellar mass, and 
then repeat the clustering measurement. The distribution 
Jj(pbar|Af*) does depend on stellar mass, but there is sub- 
stantial scatter, so it is unclear a priori what the resulting 
mark correlation function would look like. The result can be 
directly compared to the original measurement (Fig.j7^i), be- 
cause by merely shuffling the marks we are not changing the 
overall mark distribution. If the resulting mark correlation 
function is similar to the original (unshuffled) one, then this 
could be interpreted as evidence that the pbar mark correla- 
tion is due to the environmental correlation of stellar mass. 
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If the resulting mark correlation is weak but significant, then 
the environmental correlation is partly due to stellar mass, 
and if there is no mark correlation, then it is not due to 
stellar mass at all. 

We performed this test using ten bins of stellar mass (of 
width 0.15 dex), and the result yields no significant environ- 
mental correlation, as shown in Figure [9^, (open triangles in 
the lower panel). Shuffling by stellar mass appears to nearly 
completely wash out the correlation, as M Pbar is consistent 
with unity in every r v bin. Note, however, that the stellar 
mass uncertainties are comparable or larger than the bin 
widths, which could artificially weaken the correlation. 

We have also performed the same test by shuffling the 
Pbar marks as a function of g — r colour (see the distribu- 
tion in Fig. |6^i), and in this case, the pbar mark correlation 
(solid triangles in Fig. |9^i) is nearly as strong as the orig- 
inal mark correlation measurement. This suggests that the 
environmental dependence of colour partially explains that 
of Pbar- The colour-shuffled correlation function does not re- 
produce either the upturn at 400 kpc//i or the downturn at 
100 kpc/fo; however, we attribute this to the shuffling pro- 
cess. 

By taking the ratio of the marked correlation functions 
in Figure [9^,, we can make an approximate estimate of the 
fraction of the environmental dependence of pbar that is ac- 
counted for by colour and stellar mass. In particular, we use 
all of the jack-knife subsamples (not just the measurements 
in the figure) to estimate this as robustly as possible, and we 
use the mean and variance of the ratio (M r — l)/(Afpbar — 1), 
where M' is either the colour-shuffled or mass-shuffled mark 
correlation function. We use the measurements over the 
range 0.1 < r p < 2.2 Mpc//i, which encompasses the envi- 
ronmental correlations within dark matter haloes. We find 
that colour accounts for 60 ± 10 % of the environmental 
dependence of Pbar, while stellar mass accounts for only 
15±2 %. This suggests that the environmental dependence of 
colour explains the majority, but not all, of the environmen- 
tal dependence of bars. Our results are not consistent with 
Lee et al. (2012a), who claim that the environmental depen- 
dence of bars disappears at fixed colour or central velocity 
dispersion, and the disagreement may be due to their use 
of different bar classifications and environment measures; in 
addition, lenticular galaxies are excluded from their sample, 
but not from ours. 



6.1.2 Normalized mark test 

It is possible that the contribution from stellar mass is larger 
than estimated above, because the masses have larger un- 
certainties than the colours. To address this, we perform an- 
other test of the stellar mass and colour contribution to the 
bar-environment correlation, which does not involve binning 
these parameters. The purpose of this test is to remove the 
environmental dependence of stellar mass or colour (see e.g. 
Cooper et al. 2010), and consequently assess the strength of 
the residual environmental dependence of Pbar ■ 

Our procedure is as follows. For every galaxy, the pbar 
mark is normalized by the mean pbar of galaxies with that 
stellar mass (i.e., we use Pbar/(f>bar| M*) as the mark). Note 
that the mean is slightly larger than the median, which is 
plotted in Figure [6] and is similarly a smooth function of 
stellar mass, so this normalization is not sensitive to the 



mass uncertainties. Then the mark distribution is rescaled 
so that it matches the overall pbar distribution (as was done 



in Sec. 5.2 1, because consistent mark distributions are re- 



quired in order to compare mark correlations. Now the new 
mark correlation function is measured, and can be compared 
to the original one. With this test, if the mark correlation 
function were close to unity, it would mean that M* ac- 
counts for most of the environmental correlation. The same 
test is also done to assess the contribution of the g — r colour- 
environment correlation, with the analogous Pbar / (pbar \q — r) 
mark. 

The result is shown in Figure[9jD. The colour-normalized 
mark correlation function is closer to unity, and therefore 
accounts for more of the bar-environment correlation. As in 
the previous section, we can estimate the relative contribu- 
tion of colour and mass to this correlation, now using the 
ratio (Mpbar — M')/ (Mpb ar — 1)> where M' is the colour- or 
mass-normalized mark correlation function. This yields an 
estimate of 60 ± 5 % of the bar-environment correlation ac- 
counted for by colour, consistent with the mark shuffling test 
in Section [6. 1.1 1 Stellar mass now accounts for 25 ± 10 %, a 
larger contribution than estimated above, but still less sig- 
nificant than colour. We conclude that the environmental 
dependence of pbar is not primarily due to that of stellar 
mass. 

Perhaps more than stellar mass, the colour is a better 
tracer of star formation (and dust content) in disc galaxies, 
which in turn is expected to be related to the likelihood 
of the galaxies having a bar (e.g., Scannapieco et al. 2010; 
Masters et al. 2010b). Redder disc galaxies with older stellar 
populations and in more massive haloes are more likely to 
have formed a stable bar; however, mergers/interactions can 
disrupt a bar, which could explain why the original pw 
mark correlation function (circle points in Fig.|9|, unlike the 
colour-shuffled one, turns toward unity at small separations. 

6.2 Colour dependence in mock galaxy catalogues 

To add to the interpretation of the colour dependence of the 
Pbar-environment correlation in the previous section, we an- 
alyze the clustering of galaxies in a mock galaxy catalogue, 
in which we add bar likelihoods with a prescription based 
on galaxy colour. 

We use the mock catalogue of Muldrew et al. (2012), 
which was constructed by populating dark matter haloes of 
the Millennium Simulation (Springel et al. 2005) using the 
halo occupation model of Skibba & Sheth (2009). The cat- 
alogue reproduces the observed luminosity function, colour- 
magnitude distribution, and the luminosity and colour de- 
pendence of galaxy clustering in the SDSS (Skibba et al. 
2006; Skibba & Sheth 2009). Central galaxies in haloes are 
distinguished from satellite galaxies, which are distributed 
around them and are assumed to follow a Navarro, Frenk & 
White (1996) profile with the mass-concentration relation 
from Maccio, Dutton & van den Bosch (2008). 

For the purposes of this work, which is focused on disc 
galaxies, we construct a sub-catalogue from this mock, by 
selecting galaxies from the colour-magnitude distribution. 
In particular, we first select galaxies with M r — 51og(/i) < 
— 19.4. We require that the luminosity function is consistent 
with the data. (Since the catalogue was constrained with 
absolute magnitudes _K"-corrected to z = 0.1, we use the 
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Figure 9. Left: Pbar mark correlation function, shuffled as a function of g — r colour (solid triangles) and stellar mass (open triangles). 
Right: Pbar/(Pbar \q~ r ) mark (solid triangles) and Pbar/ (Pbarl-^*) mark (open triangles) correlation functions. For comparison, the original 
Pbar mark correlation function is also shown (open circles, same as Fig. [7^). Ten bins were used, most with ~1000-1400 galaxies/bin; 
using fewer or more bins yields similar results. The unweighted correlation function (upper panel in previous two figures) is omitted, 
because the full sample is used for all three mark correlation measurements. 



01 M r luminosity function.) We randomly select galaxies (in- 
dependently of halo mass or central/satellite status) in abso- 
lute magnitude bins until the consistent luminosity function 
is obtained. Secondly, we similarly require that the colour- 
magnitude distribution, p(g — r\M r ), is consistent with the 
data, using bins of 0.25 mag (see e.g. Skibba & Sheth 2009). 
The selection of disc galaxies in the GZ catalogue means 
that the red sequence is under-represented (Mil). Finally, 
we use p(pbar) (Fig. [jj), and p(pbar|# - r) (Fig. [6} distribu- 
tions to generate "pbar" for the mock galaxies. That is, we 
assume that the environmental dependence of pbar is due to 
that of g — r colour, which in turn is due to more massive 
haloes in dense environments. 

We can now measure the projected correlation func- 
tion and pbar marked correlation function of the mock cata- 
logue, in order to compare to the GZ measurements in Fig- 
ure [T] The result (averaged over eight realizations) is shown 
in Figure [To] As with the observational measurements, the 
errors are estimated using jack-knife resampling; the vari- 
ance of the eight mocks is much smaller. If we were to ap- 
ply the observed errors instead (and account for the differ- 
ent size of the GZ and mock catalogues), we obtain similar 
error bars at large scales but smaller ones at small scales 
(r p < fewMpc/ft). 

In the upper panel, the discrepancy between these pro- 
jected correlation functions at large scales has been previ- 
ously observed and is not statistically significant (see Zehavi 
et al. 2005; Skibba et al. 2006); it is likely due to cosmic 
variance. The discrepancy at small scales, however, is signif- 
icant. The fact that the correlation functions are consistent 
at scales of r v > 1 Mpc/ft, but the small-scale clustering of 
the GZ catalogue is suppressed, could mean that the satellite 
distribution as a function of halo mass is slightly different 
in the real universe, and is not reproduced with the colour- 
magnitude selection procedure. 

The Pbar(<? — t) mark correlation function of the mock 
is weaker than the GZ measurement, but similar to the 
(9 ~ reshuffled mark measurement in Figure gk. This sug- 
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Figure 10. Clustering comparison between the Galaxy Zoo and 
mock catalogues. Solid red circles indicate the projected corre- 
lation function and mark correlation function of the mock cata- 
logue, using Pbar (9 ~ r ) as the mark (see distribution in Fig.[6|i, 
and text for details). For comparison, the original GZ measure- 
ments are also shown here (open circles, same as Fig. [7ji); the 
points are slightly offset in the lower panel, for clarity. 



gests that part, but not all, of the environmental depen- 
dence of pbar is due to more massive haloes hosting redder 
galaxies, which are more likely than average to be barred. 
By taking the ratio of the marked correlation functions, 
(M moc k — 1)/(Mgz — 1), we estimate that the colour-halo 
mass correlation accounts for 50 ± 10 % of the environmen- 
tal dependence of Pbar, which is slightly lower than, but con- 
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Figure 11. Mean of halo occupation distributions in the mock 
catalogue. Black solid histogram: mean occupation number of 
all (colour-magnitude selected) disc galaxies in the mock; blue 
dashed histogram: mean occupation number of central galaxies; 
green short-dashed histogram: mean occupation number of satel- 
lite galaxies. The peak of the halo mass distribution is indicated 
by the arrow. The red dotted histogram shows the mean occu- 
pation number for galaxies with pbar > 0.2 (where the bar likeli- 
hood is computed from the p(pbarl<? — r) distribution), which are 
a combination of central galaxies in low-mass haloes and satel- 
lites in massive haloes. A larger Pb ar threshold yields a lower 
(WbarlA^halo), but of similar shape. 



sistent with the estimate in Section [6. 1| conversely, the rest 
(also 50 ± 10 %) is due to other processes unrelated to colour 
or stellar mass, perhaps involving the gas content (see Mas- 
ters et al. 2012) or angular momentum distribution. 

Also note that, as in Figure [9] the pt,ar(# — r) mark cor- 
relation function in Figure [10] lacks a drop in strength at 
r p ~ lOOkpc/ft, which we see in the original clustering mea- 
surement. This implies that, in the real universe, although 
galaxies at small separations (usually center-satellite galaxy 
pairs) tend to be redder in more massive haloes, this does not 
entail a higher bar fraction; the lack of a pbar-environment 
correlation at small separations in Figure [7|i is not related 
to galaxy colour. 

Finally, we have computed the halo mass distribu- 
tion and halo occupation distribution of galaxies in the 
mock catalogue. The halo occupation distribution (HOD) 
is the number distribution of galaxies occupying haloes of 
a given mass, and of particular importance for galaxy clus- 
tering is the mean occupation function, (N\M) (which is 
described further in Section 6.3 1. The mean occupation func- 



tions of galaxies in the mock are shown in Figure |11| The 
mock galaxies are mostly hosted by haloes with masses 
M vir > 4 x 10 11 h' 1 M©; there are fewer haloes less massive 
than this, due to the luminosity threshold (M r — —19.4). 
The central galaxy HOD drops off at high masses be- 
cause the central galaxies of these haloes rarely meet the 
CMD selection criteria of our GZ catalogue; to wit, many 



centrals in massive haloes are elliptical, not disc, galax- 
ies (Skibba et al. 2009; Guo et al. 2009; De Lucia et al. 
2011). Satellite galaxies dominate in number at masses of 
M v i r > 10 13 hT 1 Mq. In the mock, the "barred" galaxies 
(determined from the p(pbar\g — r) distribution), indicated 
by the dotted histogram, consist of a combination of cen- 
tral galaxies in low-mass haloes and satellites in massive 
haloes. The fraction of barred galaxies in the mock is not 
strongly halo mass dependent, but it is highest between 
lO 12 ' 5 < M vir < 10 14 ' 3 hT 1 Mq, in the haloes that typi- 
cally host galaxy groups. The HOD statistics of the mock 
catalogue will be compared to the results of halo occupation 
modeling, in the following section. 



6.3 Halo occupation modeling of the clustering 
measurements 

In this section, complementary to the mock catalogue anal- 
ysis of the previous section, we apply dark matter halo mod- 
els to the measured projected correlation functions, w p (r p ), 
of the whole volume-limited sample of (disc) galaxies, and 
of the subsamples of barred and unbarred galaxies, plotted 
in the upper panels of Figure [7] and Figure [SJa. Since there 
are only small differences between these measurements for 
barred and unbarred galaxies, one can expect small differ- 
ences between the well-fitting models. The purpose of the 
halo model analysis is to constrain the types of haloes that 
host barred and unbarred disc galaxies. 

We use a halo occupation model of galaxy clustering, 
(e.g., Zheng et al. 2007; Zehavi et al. 2011), in which the 
halo occupation distribution, P(N\M), of central and satel- 
lite galaxies depends on halo mass, M, and the luminos- 
ity threshold, Z/ m in- In this case the luminosity threshold 
is M r < —19.4, corresponding to an approximate halo mass 
threshold of M m in ~ 4 — 5 x 10 11 h' 1 Mq (which is consistent 



with the mock catalogues in Section 6.2 l 



Haloes of mass M are occupied by Agai galaxies, con- 
sisting of a single central galaxy and A sa t satellite galaxies, 
such that the mean occupation function is described as the 
following: 

(A ga i|M,L min ) = <A con jM,L min ) [l + <A S at|M,L min >] (8) 
where, 



(Ncen\M) = - 



1 + erf 



/ log(M/M min ) 

\ OlogM 



and 



(JV aat |Af> = 



M - M 
M[ 



(9) 



(10) 



In practice, we account for the fact that there is significant 
scatter in the relation between central galaxy luminosity and 
halo mass, and that the satellite halo occupation function 
drops off more rapidly than a power-law at low masses just 
above M min . See Appendix A2 of Skibba & Sheth (2009) for 
details. 

We will also use the halo occupation models to compare 
the fraction of satellite galaxies of barred and of unbarred 
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Figure 12. Halo occupation models of clustering of all galaxies in 
the sample (solid black curve), the barred subsample of galaxies 
(red dashed curve), and the unbarred subsample (blue dotted 
curve). The corresponding measured correlation functions are also 
shown, for barred and unbarred galaxies (solid and open triangle, 
respectively) . 



sample 



log M Inia /h- 1 M Q M 1 /M u 



all 

barred 
unbarred 



11.65 ±0.05 
11.65 ± 0.05 
11.65 ± 0.05 



16.0 ± 1.0 0.90 ±0.04 
16.0 ± 1.0 0.95 ±0.06 
19.0 ± 1.0 0.80 ±0.06 



Table 1. Halo occupation distribution parameters for the model 
| |8)10| fitted to the clustering measurements of the whole sample 
and the barred and unbarred subsamples, shown in Fig. 1 1 2| 



galaxies. The satellite fraction is given by 

/ . dM (dn/dM) (N sat \M) 



./"sat 



f M . dM (dn/dM) «A ce „|M) + <JV sat |Af» ! 



(11) 



where dn/dM is the halo mass function (Sheth & Tormen 
1999; Tinker et al. 2008b). Note that we will not attempt to 
account for the fact that central galaxies in massive haloes 
will often not meet the selection criteria for disc galaxies, 
because these galaxies will be dominated in number by satel- 
lites (see Fig. whose abundance we can constrain. 

In Figure [12] we show the results of the halo occupa- 
tion modeling, applied to the whole catalogue and to the 
subsamples of barred and unbarred galaxies. The parame- 
ters of these HOD models are listed in Table [T] We fixed the 
parameters <Ti og Af = 0.26 and Mo = 10 11 ' 60 Mq/Ii, because 
they are not constrained well by HOD models of cluster- 
ing (see Zheng et al. 2007; Zehavi et al. 2011); (TiogM, which 
quantifies the scatter between central galaxy luminosity and 
halo mass, can be constrained by satellite kinematics and the 
conditional luminosity function (More et al. 2009; Cacciato 
et al. 2009). 

As stated above, the halo mass threshold of the three 



measurements is the same. Nonetheless, because of differ- 
ences in the small-scale clustering, there are differences in 
the satellite HOD, (7V sat |M) (Eqn.[l0|. In particular, firstly, 



the fraction of satellite galaxies varies. The whole catalogue 
has / sat ~ 23%, consistent with the mock catalogue analysis 
in Section 6.2 which yielded a similar fraction (also ~ 23%). 



For comparison, the barred and unbarred subsamples have 
/ sa t « 25% and 20%, respectively. 

Secondly, the key difference between the well-fitting 
models for the barred subsample is that they have a steeper 
slope a (compared to the slope for the full sample and for 
unbarred galaxies), which means that the larger satellite 
fraction of barred galaxies is due to more satellites in more 
massive haloes. In contrast, the well-fitting models for the 
unbarred subsample have a shallower slope, so not only is 
the unbarred sample dominated by central galaxies in lower 
mass haloes, but the small fraction of unbarred satellites is 
not in the most massive haloes either. 



7 CONCLUSIONS AND DISCUSSION 

We selected a volume-limited catalogue of 15810 nearly face- 
on disc galaxies in the SDSS, which have visual morphology 
classifications from Galaxy Zoo 2. We analyzed the proper- 
ties of galaxies with bars and bulges, characterizing bar and 
bulge likelihood with the pbar and fracdeV parameters. Us- 
ing "marked" two-point correlation functions, we quantified 
the environmental dependence of bar and bulge likelihood 
as a function of the projected separation between galaxies. 

To conclude, the following are the main results of our 
paper: 

• Correlations of bars and bulges with colour and 
stellar mass: We find a strong correlation between the bar 
likelihood (pbar) and optical colour and stellar mass, such 
that redder and more massive disc galaxies are up to twice 
as likely to have bars than their bluer low-mass counterparts, 
although there is considerable scatter in the correlation, es- 
pecially at the red (high-mass) end. We find similar corre- 
lations with bulge strength (fracdeV), but with less scatter. 
The quantities pbar and fracdeV appear to have a transition 
at the same stellar mass and colour (M* « 2 x 10 10 M Q , 
g - r w 0.6). 

• Environmental dependence of bars and bulges: 

We clearly detect and quantify the environmental depen- 
dence of barred galaxies and of bulge-dominated galaxies, 
such that barred and bulge-dominated disc galaxies tend 
to be found in denser environments than their unbarred 
and disc-dominated counterparts. In particular, by analyz- 
ing pbar and fracdeV marked correlation functions, we ob- 
tained environmental correlations that are statistically sig- 
nificant (at a level of > 6a) on scales of 150 kpc to a few 
Mpc. From sparse sampling tests with our catalogue, we ar- 
gue that the small number statistics of previous studies in- 
hibited their detection of a bar-environment correlation. 

• Contribution from colour and stellar mass to 
bar-environment correlation: By accounting for the en- 
vironmental dependence of colour and stellar mass, we argue 
that they contribute approximately 60±5 % and 25±10 %, of 
the pbar-environment correlation, respectively. From a sim- 
ilar analysis of a mock galaxy catalogue, we argue that the 
environmental dependence of Pbar appears to be partially 
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(50 ± 10 %) due to the fact that redder galaxies, which are 
often barred, tend to be hosted by more massive haloes. 
Conversely, up to half of the bar-environment correlation is 
not due to colour or stellar mass, and must be due to envi- 
ronmental influences or to another independent parameter 
(possibly gas content, or angular momentum distribution). 

• Halo model analysis of clustering of barred 
galaxies: Our analyses with a mock galaxy catalogue and 
halo occupation models suggest that barred galaxies are of- 
ten either central galaxies in low-mass dark matter haloes 
(M ~ 10 12 Mq) or satellite galaxies in more massive haloes 
(M ~ 1O 13_14 M , hosting galaxy groups). 

We argue that the environmental dependence of galaxy 
colours can account for approximately a half (50 ± 20 %) 
of the environmental correlation of Pbar- The optical colours 
are correlated with star formation rate and age, as well as 
gas and dust content, all of which may be related to the 
presence of disc instabilities such as bars. We find that a 
galaxy's stellar mass and bulge component, on the other 
hand, do not appear to be strongly related to its likelihood 
of having a bar. This suggests that it is primarily older disc 
galaxies with lower star formation rates (which often reside 
in denser environments) that are able to form and maintain 
a stellar bar (see Masters et al. 2010b). 

Conversely, this means that the remaining half of the 
environmental correlation o/pbar is not explained by the en- 
vironmental dependence of colour or stellar mass, suggesting 
that bar formation (or the lack of bar destruction) is likely 
influenced by the galaxy's environment, in addition to the 
effects described above. Bulge formation, which is to some 
extent independently correlated with the environment (see 



Section 5.2|, is also expected to be affected by interactions 



and merger activity (Hopkins et al. 2009; Kannan et al., in 
prep.). 

During the final stages of this work, Martinez & Muriel 
(2011) in a related study found that the bar fraction does 
not significantly depend on group mass or luminosity, or 
on the distance to the nearest neighbour. Their sample is 
smaller than ours, however, and is apparent magnitude- 
limited rather than volume-limited. In addition, they use 
bar classifications from Nair & Abraham (2010), which in- 
clude somewhat weaker bars than Galaxy Zoo 2 (see Mil), 
which are bars that tend to be found in bluer galaxies (and 
hence in less dense environments). In another recent paper, 
Lee et al. (2012a) also analyze the environmental depen- 
dence of bars, using bar classifications consistent with Nair 
& Abraham (2010), and claim that the bar fraction does not 
depend on the environment at fixed colour or central veloc- 
ity dispersion, contrary to our results. However, a crucial 
difference between these two analyses and ours is that they 
use environment measures that mix environments at differ- 
ent scales, while we analyze the environmental correlations 
as a function of galaxy separation. 

A particularly interesting result of this paper is the 
scale dependence of the environmental correlations of bar 
and bulge likelihood (see Fig.[7j Sec.[5|. Environmental cor- 
relations should be interpreted differently at different scales 
(Blanton & Berlind 2007; Wilman et al. 2010; Muldrew et al. 
2012), as galaxies at small separations (r p < 2 Mpc//i) are 
often hosted by the same dark matter halo, while galaxies at 
larger separations are often hosted by separate haloes. We 



see that more massive haloes, which tend to reside in rela- 
tively dense environments, tend to host more disc galaxies 
with bars and bulges. 

Moreover, the pbar-environment correlation peaks at 
r v ~ 400 kpc//i, which suggests that many barred galax- 
ies are central or satellite galaxies in groups and clus- 
ters. That is, some aspect of group environment triggers 
the formation of bars, in spite of the fact that bars are of- 
ten thought to form by internally driven secular processes 
(e.g., Kormendy & Kennicutt 2004). Secular processes may 
sometimes be externally driven. For example, cosmologi- 
cal simulations predict that tidal interactions with dark 
matter substructures, which are common in such environ- 
ments, could induce bar formation and growth (Romano- 
Diaz et al. 2008; Kazantzidis et al. 2008). On the other hand, 
the pbar-environment correlation is not significant for closer 
pairs, suggesting that the enhanced likelihood of galaxies 
being barred is erased if the galaxies are merging with 
each other; however, this measurement at small separations 
(r p < 100 kpc/h,) has large uncertainty and may be affected 



by fiber collisions (see Sec. 2.4 1, so it should be viewed with 



caution. Analyses of bars in close pairs of galaxies (e.g., Nair 
& Ellison, in prep.) could shed more light on this issue. 

In general, we can at least conclude that group environ- 
ments increase the likelihood of bar formation in disc galax- 
ies. Minor mergers and interactions are relatively common 
in galaxy groups (Hopkins et al. 2010b), and tidal interac- 
tions with neighboring galaxies can trigger disc instabilities 
and subsequent bar formation (Noguchi 1996; Berentzen et 
al. 2004); there is observational evidence for this as well 
(Elmegreen et al. 1990; Keel et al. 1996). Tidal interactions 
can also affect the bar's pattern speed and other properties 
(Miwa & Noguchi 1998). The evolution of barred galaxies 
in group environments and in minor mergers/interactions is 
clearly in need of further study. 

Considering that bar formation does appear to depend 
on the host galaxy's environment, and that bars form by 
secular evolution, our results suggest that the dichotomy be- 
tween internal secular processes and external environmental 
processes is not as strict as previously thought. It is possible 
that some structural changes in galaxy discs may be trig- 
gered or influenced by the galaxy's environment. For exam- 
ple, Kormendy & Bender (2012) recently argued that, "ha- 
rassment", the cumulative effect of encounters with satel- 
lite galaxies, may influence secular evolution. Furthermore, 
"strangulation", in which the hot diffuse gas around newly 
accreted satellites is stripped, removes the fuel for future 
star formation (Larson et al. 1980), and could contribute to 
more stable or growing bars (Berentzen et al. 2007; Masters 
et al. 2012). 

In addition, our results could also indicate a link be- 
tween bars and active galactic nuclei (AGN) activity. We 
have shown that barred galaxies tend to reside in dense 
group environments, while galaxies hosting AGN also tend 
to be found in such environments (e.g., Mandelbaum et al. 
2009; Pasquali et al. 2009). In some models, AGN are as- 
sumed to be fueled by recent mergers; however, some have 
also argued that bars and disc instabilities may be an inter- 
nal mechanism through which low angular momentum gas 
is driven towards the nucleus (Bower et al. 2006; Hopkins & 
Quataert 2010; McKernan et al. 2010). Nonetheless, a cor- 
relation between barred galaxies and AGN activity has not 
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been detected observationally (Lee et al. 2012b; Cardamone 
et al., in prep.). 

Lastly, we note that our results can be used to constrain 
galaxy formation models, such as the semi-analytic models 
of Benson & Devereux (2010) and De Lucia et al. (2011), 
and the hydrodynamic simulations of Heller et al. (2007), 
Croft et al. (2009), Sales et al. (2012). Marked correlation 
functions, and marked clustering statistics in general, are 
sensitive to environmental correlations at different scales, 
such that small changes in model parameters could yield 
environmental dependencies of galaxy bars and bulges that 
can be compared to measurements with Galaxy Zoo (see 
Figures [7]|9|. In addition, a result from our halo model anal- 
ysis is that barred galaxies tend to be central galaxies in 
lower mass haloes (Mhaio ~ 10 12 Mq) and satellite galaxies 
in more massive haloes (Mhalo ~ 10 13 Mq), which can also 
be compared to other models. 
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APPENDIX A: ERRORS OF THE 
CLUSTERING MEASUREMENTS 

As discussed in Section[3] we use jack-knife resampling to es- 
timate the statistical errors of our clustering measurements 
w p (r p ) and M p (r p ) (Figures 7j9l. The jack-knife errors of 
our marked projected correlation functions tend to be larger 
than the Poisson errors (not shown). We note that Norberg 
et al. (2009) has shown that the jack-knife method does not 
recover the scale dependence of errors of the (unmarked) 
correlation function, often overestimating the errors at small 
scales, and the results are sensitive to the number of sub- 
catalogues into which the data is split. Although our un- 
certainty estimates are important, our primary results are 
not particularly sensitive to the precise value of the errors. 
We have performed additional error analyses with twice as 
many jack-knife sub-catalogues, and obtained similar errors 
(within 10%) at all scales for both the correlation functions 
and marked correlation functions. 

We estimate the error in each r p (projected galaxy sep- 
aration) bin by computing the variance of the measurements 
of the jack-knife sub-catalogues. For example, clustering 
measurements of the thirty jack-knife sub-catalogues used 
to estimate the errors of the measurements in Figure [7^, are 
shown in Figure |A1| Most of the jack-knife sub-catalogues 
yield similar measurements, although a single outlying mea- 
surement is responsible for a significant fraction of the er- 
ror. Outliers may be due to anomolously large structures or 
voids, or in the case of mark clustering, to noise in the marks. 
Nonetheless, the outlying measurement is not an extreme 
outlier, and the Sloan Great Wall, an extremely massive 
superstructure that can influence clustering measurements 
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Figure Al. Jack-knife clustering measurements for Fig. [7^.. 

(Zehavi et al. 2011; Norberg et al. 2011) is beyond our up- 
per redshift limit (z < 0.060). The jack-knife measurements 
for the lower panel of Figure [TJa (with the fracdeV mark) 
are similar, also with a single outlier. If the outlier were 
excluded, the resulting jack-knife error estimates would be 
lower by 21% on average for w p , by 18% for Mpbar, and by 
25% for Mfracdcv- 

We also compute jack-knife covariance matrices 
(Eqn. [6| for each of the clustering measurements. We show 
the covariance matrix of the measurement in the lower panel 
of Figure [7^, (pbar mark) in Figure |A2| Most of the errors 
are not strongly correlated, although the r p bins centered at 
1 and 2 Mpc/h are weakly correlated with smaller-scale bins 
(most of which have stronger mark correlations) , which may 
explain why the mark correlation measurements of these 
bins are slightly larger than in the neighboring r p bins. 

In any case, the range of jack-knife clustering measure- 
ments in Figure |A1| is not extremely large, and the corre- 
lations between the r p bins in Figure [A2| are not extremely 
strong, so we conclude that the measured mark correlation 
functions are robust, as are the inferred environmental cor- 
relations. 
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10 Mpc/h 
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Figure A2. Jack-knife covariance matrix for the Mpt, a r( r p) mea- 
surement of Fig. [7p,. The grayscale as a function of galaxy sep- 
aration indicates regions at which the measurement errors are 
correlated, such that darker regions have more correlated mea- 
surements. 
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